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ABSTRACT 

The nature of the magnetic coupling between T Tauri stars and their disks determines 
not only the mass accretion process but possibly the spin evolution of the central star. 
We have taken a recently-published surface magnetogram of one moderately-accreting 
T Tauri star (V2129 Oph) and used it to extrapolate the geometry of its large-scale 
field. We determine the structure of the open (wind-bearing) field lines, the closed (X- 
ray bright) field lines and the relatively small subset of field lines that pass through 
the equatorial plane inside the Keplerian co-rotation radius and which are therefore 
available to accrete. 

We consider a series of models in which the stellar magnetic field is opened up 
by the outward pressure of the hot coronal gas at a range of radii or source surfaces. 
As the source surface is increased, accretion takes place along progressively simpler 
field structures and impacts on progressively fewer sites at the stellar surface. This is 
consistent with the observed variation in the Ca II IRT and Hel lines which suggests 
that accretion in the visible hemisphere is confined to a single high-latitude spot. By 
determining the density and velocity of the accretion flows, we find that in order to 
have most of the total mass accretion rate impacting on a single high-latitude region 
we need disk material to accrete from approximately 7R*, close to the Keplerian co- 
rotation radius at 6.8R*. 

We also calculate the coronal density and X-ray emission measure. We find that 
both the magnitude and rotational modulation of the emission measure increase as the 
source surface is increased. For the field structure of V2129 Oph which is dominantly 
octupolar, the emission forms a bright, high-latitude ring that is always in view as 
the star rotates. Since the accretion funnels are not dense enough to cause significant 
scattering of coronal X-ray photons, they provide only a low rotational modulation of 
around 10% at most. 
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1 INTRODUCTION 

T Tauri stars are young solar-like stars that are still in the 
process of contracting onto the main sequence. The so-called 
classical T Tauri stars (cTTS) show signs of active accretion 
and are typically younger than their counterparts the weak- 
line T Tauri stars. One of the enduring puzzles about cTTS 
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is their slow to moderate rotation, which suggests that they 
must have succeeded in losing angular momentum as they 
contracted. This can in principle be achieved by magneti- 
cally coupling the star to its disk and indeed the presence 
of a disk, as indicated by IR excesses does seem to be corre- 



lated with slower rotation (e.g. Rebull et al. (20061). While 



the detailed mechanisms by which this coupling is achieved 
are still unclear, the general hypothesis is well established. 
The stellar magnetic field disrupts the inner edge of the disk, 
forcing accreting material to flow along field lines down to 
the stellar surface. The torque associated with this magnetic 
field may balance the accretion torque, thus preventing the 
star from spinning up. Any attempt by the star to spin up 
generates a stronger magnetic field which exerts a greater 
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Figure 1. The variation with radius of (solid) the radial compo- 
nent of the magnetic field and (dotted) the non-radial component. 
In both cases, the value plotted is averaged over latitude and lon- 
gitude. Lines are drawn for a range of source surfaces. At the 
source surface itself, both Bg and Ba, are identically zero. 



torque that acts to slow the star down. Field strengths at 
the stellar surface can then be estimated by calculating the 
values at the disk necessary for such a torque balance and 
extrapolating back to the surface, assuming some simple ge- 



ometry (Konigl 1991 Collier Cameron & Campbell 1993 



Shu et al. 1994 1 . This process, however, requires surface field 
strengths of around 0.1-5 kG and a disk that is truncated 
close to the Keplerian corotation radius. 

The models can in principle be tested by measuring the 
surface field strength of T Tauri stars. Zeeman-Broadening 
measurements can now be used to measure mean fields 



on the stellar surface (Johns-Krull et al. 2004| Valenti & 
| Johns-Krull|2004| | Yang et al.|2005||Johns-Krull|2007| ). This" 

gives distributions of local field strengths on the stellar sur- 
face that may be as high as 6kG which is well within the 
range predicted by theory. However, these are similar on 
both accreting and non-accreting T Tauri stars, suggesting 
that disk-locking is not responsible for determining the field 
strength. There also seems to be little correlation with either 
the values predicted by theory, or with the stellar rotation 
rotation rate or Rossby number. Furthermore, although the 
average surface field strengths reach several kG, current ob- 
servations indicate that the strength of the dipole compo- 
nent is at least an order of magnitude below what is required 
for disc-locking models to operate (e.g. Daou, Johns-Krull 



fc Valenti|2006| |Yang, Johns-Krull fc Valenti|20"07 1 

A complementary approach is to use circular polari- 
sation measurements both in photospheric lines (presum- 
ably formed across the whole stellar surface) and the He- 
lium A5876A emission line believed to form in the accretion 
shock ( jBeristain et al.||2001| |Yang et aL|2007| ). Polarisation 
measurements taken at a range of rotation phases can then 
be used to provide information on the structure of the mag- 
netic field both across the surface and in the accretion shock. 
The low circular polarisation measured in photospheric lines 
rules out the presence of a global dipole, suggesting per- 
haps a complex field that may be locally very intense, but 
which is organised into small enough scales that the Zee- 
man signatures are complex in shape and associated with 
small longitudinal fields. In comparison, the large, rotation- 
ally modulated unipolar Zeeman signatures in the Helium 




Figure 2. The coronal magnetic structure of V2129 Oph for a 
source surface at 3.2R* (top) and at 6.8R* (bottom). Closed field 
lines are shown white, open field lines are shown blue. 



emission lines suggest that the field associated with accre- 
tion is much simpler. This is no doubt to be expected, since 
it is the largest-scale field lines that will interact with the 
disk and these are likely to have the simplest structure. The 
nature of the star-disk interaction will clearly depend on the 
structure of the whole coronal field, and in particular the 
degree to which the complexity that seems apparent at the 
stellar surface may have died away by the time the coronal 
field reaches the inner edge of the disk. 

Signatures of this complexity may be apparent in the 
X-ray emission from T Tauri stars. In addition to the ubiqui- 
tous flaring-like variability seen in the X-ray emission of the 
Chandra Orion Ultradeep Project (COUP; Getman et al. 
2005 ) stars, there is also a clear rotational modulation in a 



small but significant number of these stars ( Flaccomio et al. 
|2005p . Such modulation has previously been difficult to de- 
tect because it requires continuous observation over several 
rotation periods ( Giidel et al.| 1995 Ma rino et al.|200"3 Hus- 
|sain et aL]|2005[ ). This rotational modulation suggests that 
the X-ray coronae of T Tauri stars are compact, with dis- 
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crete bright and dark regions. Some contribution to the rota- 
tional modulation may also come from the accretion funnels 
themselves, which, if they are dense enough, may obscure 
the underlying emission. Gregory et al. ( |2006[ ) have mod- 
elled the structure of T Tauri coronae and found significant 
rotational modulation, with X-ray periods that are typically 
equal to the optical period (or one half of it). While there 
seems to be little time-correlation between optical and X- 
ray variability, there appears to be a significant correlation 
between optical variability and X-ray luminosity. This sug- 
gests that the X-ray emission of T Tauri stars is primarily 
produced in a complex arrangement of magnetically-active 



regions in a manner similar to observed on the Sun ( Stassun 
et al.|2006||2007 l 



One singularly non-solar aspect of T Tauri coronae how- 
ever is the possible interaction between the magnetic field 
of the star and the accretion disk. Both in the COUP study 
and also more recently in the DROXO (Deep Rho Ophi- 
uchi Xmm-Newton Observation) study, there are flares for 
which modelling suggests the loops responsible must be large 



perhaps reaching out into the accretion disk ( Favata et al 



20051. Field lines that connect the star and the disk in 



this way will be stressed by their interaction with the disk 



and may indeed carry much of the accretion torque (von 
Rekowski fc BrandenbmjpXM] |Long et al.|[2007| . The na- 
ture of the star-disk interaction may therefore depend on 
the structure of the magnetic field that links to the disk, 
whether it is primarily open as in these more recent models, 



or closed as in the traditional disk-locking models (Konigl 
199lj |Collier Cameron fc Campbell||l993l |Shu et al.||1994| 



Romanova et al.|2002||2003||Long et al.|2005| >. Recent mod- 
els of accretion powered stellar winds suggest that the an- 
gular momentum gained by the star due to accretion may 
be lost in a stellar wind, rather than exported back into the 
disk ( |Matt fc Pudritz|2008a|b"l ). The degree of complexity of 
the field affects the torque that it can exert and the location 
of the inner edge of the accretion disk. 



Recent results by Jardine et al. ( 2006 1 suggest that 



the nature of this star-disk interaction may be fundamen- 
tally different in high and low mass stars. If we consider 
pre-main sequence stars of progressively lower mass, their 
pressure scale heights A will increase (since for a polytrope 
R oc A4" _1//3 and so A k M* 5 ^ 3 ) while their co-rotation radii 
will decrease. Thus the lower the stellar mass, the more likely 
it is that the co-rotation radius will be within the corona. 
If a disk is present, this may mean that the disk strips off 
the outer edge of the corona, converting field lines that were 
closed and therefore possibly X-ray bright, to an open, X- 
ray dark, wind-bearing configuration. Even in the absence 
of a disk, of course, centrifugal effects may do the same job, 
preventing the corona from extending much beyond the co- 
rotation radius (Saficr 1998). This may well explain the drop 
in X-ray emission measures towards lower masses seen in the 
COUP dataset (Getman et al. 2005) but it also suggests that 
for the higher-mass stars, whose coronae may not extend as 
far as the co-rotation radius, any disk that exists at the 
co-rotation radius may interact not with the closed X-ray 
emitting field lines of the star but rather with the open field, 
which in the absence of a disc would carry a stellar wind. It 
is not clear what effect this would have on the interaction of 
the stellar and disk magnetic fields, and in particular on the 
location of the inner edge of the accretion disk. This may be 



determined by a balance between the matter energy density 
of the disk and the magnetic energy density of the stellar 
magnetic field, or by a small-scale balance between the disk 
viscosity and the resistivity which allows matter to transfer 



from disk field lines to stellar field lines (e.g.|Romanova et al 



(2002 1; Collier Cameron & Campbell ( 1993|). It seems clear 
that not only the strength of the stellar field but also its 
structure are important. Recent models by Bcssolaz et al. 
(2008) suggest that accretion funnels may indeed form in 
mildly-accreting T Tauri stars with weak dipolar fields of 
only 140G, but that higher field strengths or different field 
topologies are needed to explain more typical accretion rates 
of M ~ 10" 8 Mgyr -1 . 

In this paper we use recently-published surface magne- 
tograms of a fairly typical T Tauri star to determine the 
structure of its corona and in particular to investigate the 
relationship between the closed-field regions that may con- 
tribute to the X-ray emission, the open field lines that may 
contribute to the wind and the field lines that pass through 
the equatorial plane inside of co-rotation which are available 
to accrete. By solving for the flow along these accreting field 
lines we determine the influence of the accretion funnels on 
the predicted magnitude and rotational modulation of the 
X-ray emission and also the mass accretion rate. 



2 THE CTTS V2129 OPH 

V2129 Oph (SR 9, AS 207A, GY 319, ROX 29, HBC 264) is 
a fairly typical, moderately accreting classical T Tauri star 



with rotation rate of 6.53 d (Shevchenko & Herbst 



Many of its parameters are derived in Donati et al 



1998) 



(20071 



and we simply summarise their results here. They find a ra- 
dius and mass of R t = 2.4±0.3 R and M 4 = 1.35±0.15 M 
(which imply that the co-rotation radius is at 6.8i?*) and us- 
ing the evolutionary models of Siess et al. ( 2000 1, they obtain 



an age of about 2 Myr and conclude that V2129 Oph is no 
longer fully convective, with a small radiative core of mass 
~ 0.1 M* and radius ~ 0.2 _R + . The presence of a radiative 
core is relevant to the nature of magnetic field generation 
since it provides a region of high shear at the boundary be- 
tween the core and the convective region within which field 
can be generated. The mass accretion rate may be derived 
from the equivalent width of t he 8662A Ca H line emission 
g iving 4 x 10" 9 M Q yr" 1 ( |Mohanty et al.|2005| ). |Eisner 
et al.| ( 2005 I suggests that M is si gnificantly larger, about 



3x 10 _B M© yr . We follow |Donati et al.| ( |2007[ ) and assume 
here that M ~ 10 -8 M s yr _i , i.e., a conservative compro- 
mise between both estimates and a typical value for a mildly 
accreting cTTSs. 

The surface magnetic field maps derived by |Donati et al.| 
( |2007[ ) are based on circular polarisation Zeeman signatures 
detected in both photospheric lines and also in emission lines 
tracing magnetospheric accretion. The surface field that is 
recovered is complex, with many mixed-polarity regions. It is 
dominated by a 1.2kG octupole, tilted by about 20° to the 
rotation axis. A weaker dipole component is also present, 
with polar strength of 0.35kG and a similar tilt of about 
30°. For the hemisphere in view, accretion appears to be 
concentrated in a single high-latitude region which covers 
about 5% of the stellar surface and is in view at around ro- 
tation phase 0.8. This apparently coincides with a dark po- 
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lar feature at the photo-spheric level. Observational evidence 
that the visible accretion spot is located at high latitudes is 
strong; the low amplitude fluctuations in the radial velocity 
of the He I line (or the Ca If IRT emission core) on the one 
hand, and the fact that the accretion spot (and associated 
Zeeman signature) is always visible to the observer on the 
other hand, both argue in favour of this interpretation. This 
implies that any accretion from the disk must flow up out 
of the equatorial plane in order to impact the stellar surface 
at high latitudes. 
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3 EXTRAPOLATING FROM THE SURFACE 
FIELD 

Once the surface magnetic field is known it is possible to ex- 
trapolate the coronal magnetic field using a variety of meth- 
ods and assumptions. Two of the most commonly-used are 
a full MHD solution and the "Potential Field Source Sur- 
face" method which was originally developed by Alt schuler] 
|fc Newkirk, Jr. | ( |1969| ) for extrapolating the Sun's coronal 
field from solar magnetograms. The latter has the advan- 
tage of computational speed and simplicity, but it cannot 
reproduce time-dependent or non-potential effects, both of 
which are likely to be important for the structure of the 
field in the region of the accretion disk. The former method 
can treat both of these effects, but is computationally very 
expensive and requires some assumption of an equation of 
state (often a polytrope is used) and the nature of the en- 
ergetics. A review and comparison of the two techniques is 
given in 



Riley et al. ( 2006 1 



While a full MHD solution is desirable in order to deter- 
mine the nature of the interaction of the stellar field with the 
accretion disk, we postpone this approach for the moment, 
and concentrate on the first question that should be asked: 
is this surface magnetogram consistent with the inferred lo- 
cation of the accretion flows? These flows originate in the 
accretion disk and if they are magnetically-channelled, the 
location at which they arrive at the stellar surface must de- 
pend on the structure of the coronal magnetic field. Only 
some fraction of the field lines that leave the stellar sur- 
face ever come close enough to the accretion disk to interact 
with it. This interaction is likely to be time-dependent and 
to inject stresses into the magnetic field, just as the con- 
vective motions, surface differential rotation and meridional 
flows stress the magnetic field lines that emerge through the 
stellar surface. Nonetheless, the subset of stellar field lines 
that pass through the disk and so are available to accrete 
is likely to be determined principally by the nature of the 
stellar field. We concentrate therefore on determining this 
overall structure and on determining which field lines are 
available for accretion and where they connect to the stellar 
surface. 

We therefore use the "Potential Field Source Surface" 
method with a code originally developed by van Ballcgooijcn 



Jardine 



et al. ( f998). Since the method has been described in^ 

et al.| ( |2002a[ ) we provide only an outline here. Briefly, we 
write the magnetic field B in terms of a flux function ^ such 
that B — — V\& and the condition that the field is poten- 
tial (V x B = 0) is satisfied automatically. The condition 
that the field is divergence-free then reduces to Laplace's 
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Figure 3. Magnetic field strength of a range of field lines at the 
point where disk material latches onto the field. If the source 
surface is beyond the co-rotation radius, this location is taken to 
be the co-rotation radius, if, on the other hand, the source surface 
is inside co-rotation, this location is taken to be the source surface 
itself. Results are shown for a range of source surface locations. 



equation V $ = with solution in spherical co-ordinates 
(r,M) 

jv i 

* = J2 [a ' mrl + b ^- {l+1) \P lm {6y m ^ (1) 

l — 1 m — — l 

where the associated Legendre functions are denoted by P; m . 
This then gives 
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The coefficients aim and bi m are determined by imposing the 
radial field at the surface from the Zeeman-Doppler maps 
and by assuming that at some height r = R s above the 
surface (known as the source surface) the field becomes ra- 
dial and hence B e (R s ) — B ( p(R s ) — 0. This second condi- 
tion models the effect of the plasma pressure in the corona 
pulling open field lines to form a stellar wind. It imposes a 
relation between the term in |l| that grows with radius and 
has coefficient ai, n and the term that decays with radius and 
has coefficient b; m , since it implies that 



aim 

bim 



Hence as the source surface radius decreases, the ratio 
aim/bim increases and therefore the contribution to ^ of the 
part that grows with r becomes relatively larger. As shown in 
Fig. |TJ, this means that at any given height, as the source 
surface is decreased, the average radial component of the 
field will increase while the non-radial contribution will de- 
crease. Changing the source surface therefore modifies the 
field structure of the whole corona, such that, as shown in 



Fig. |2j, increasing the source surface lead to field lines that 
have a smaller radial component and are therefore rounder. 

For each given surface map, therefore, the source sur- 
face is the one free parameter in determining the structure 
of the coronal magnetic field. Whatever value is chosen for 
this, however, the dominantly octupolar nature of the mag- 
netic field close to the surface remains, as does the transition 
to a dipolar field at greater heights. This transition can be 
seen in Fig. [T] where the radial field falls off with height as 
r~ 5 (appropriate for an octupolar field) close to the star 
around 2i?*, while a dependence closer to r~ 3 (appropriate 
for a dipole) becomes apparent only at heights greater than 
around Given the strengths of the dipolar and octupo- 
lar terms, we would expect the dipole to dominate at heights 
above 2R+. 
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4 THE ACCRETING FIELD 

Once we have determined the coronal structure of the mag- 
netic field, we can divide the field lines into three types, 
depending on whether they are open, closed or possibly ac- 
creting. The only field lines that have the opportunity to 
accrete are those that pass through (or close to) the disk 
at a radius where the gravitational acceleration is directed 
towards the star. Field lines that never pass near to the 
disk region may either be open (i.e. wind-bearing) or closed 
(i.e. X-ray bright). Clearly, with our field extrapolation tech- 
nique we cannot determine the fate of the field lines that 
pass through the disk. The nature of their interaction with 
the disk material will depend not only on the global disk 
properties, such as the mass accretion rate, but also on the 
small-scale physics of the plasma diffusivity and viscosity. 
An important factor is also the strength and geometry of 
the field lines that intersect the disk. We show in Fig. (|3| 
the range of field strengths of those field lines that pass 
through the disk. Models with a range of source surfaces are 
shown, some well within the Keplerian co-rotation radius 
and some beyond it. For the models where the source sur- 
face is inside the co-rotation radius, material accretes from 
the source surface along open field lines. In this case, the 
fall-off in field strength is just what would be expected from 
Fig. Q. For models where the source surface is beyond the 
co-rotation radius, however, we note that material can only 
accrete from within the co-rotation radius. In this case, ac- 
cretion takes place from a fixed radius along predominantly 
closed field lines whose field strength now varies much less 
with the value of the source surface. Field lines that pen- 
etrate the disk beyond the co-rotation radius are likely to 
be sheared open, adding to the flux of field lines that have 
been opened up by the pressure of the hot coronal gas. In 
this case our model would underestimate the flux of open 
field. 

Considerable progress has already been made in study- 
ing the interaction of simulated stellar magnetic fields with 
accretion disks ( von Rekowski &: Piskunov|2006 Long et al 



2007||Matt fc Pudrit z 2008a b). Recent models of multipolar 
fields in particular show that accreting material may take 
the most energetically favourable path (i.e. closest to the 
equatorial plane) and so the form of the accretion funnels 
depends sensitively on the stellar field structure | |Romanova| 
et al.|2 007). The distortion of the field that will take place as 
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Figure 4. Density profiles (solid lines) along example field lines. 
The source surface is at 2.5R* (top) or 6.8R* (bottom). The 
dashed lines show, for comparison, density profiles along a dipole 
field. 



it couples to the disk material is, however, beyond the scope 
of this paper and will be the subject of future work. Our 
present aim is to determine as a starting point which field 
lines are capable (by virtue of their location) of connecting 
to the disk at all. 

We therefore consider field lines that pass through the 
disc within the co-rotation radius. We assume that the disc 
lies in the equatorial plane and is of opening angle 10°. The 
interaction of the stellar magnetic field with the disc may 



lead to a warping of the inner disc ( Bouvier et al. 2006 1 



in which case the disc may lie along the magnetic equator, 
rather than the stellar rotational equator. In fact, we have 
found that allowing accretion to start from the magnetic 
equator rather than the rotational equator makes little dif- 
ference in this case, and so in the absence of any evidence 
for a warped disc in V2129 Oph, we consider the disc to lie 
in the plane defined by the stellar rotation axis. Along the 
field lines that pass through this plane inside the co-rotation 
radius, we determine the density and velocity structure of 
the accretion flow. 



Following Gregory et al. ( 2007 1 we assume that the disc 



supplies mass at a constant rate. We use a spherical grid and 
assume that each grid cell within the disc which is accreting 
supplies a mass accretion rate that is proportional to its 
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Figure 5. Average infall speed (solid line) as a function of source 
surface radius. The dashed lines show the upper and lower limits 
of the range of infall speeds. 



surface area. For example, if an accreting grid cell has a 
surface area that is 1% of the total area of all accreting grid 
cells, then this grid cell is assumed to carry 1% of the total 
mass accretion rate supplied by the disc. Therefore, if grid 
cells which constitute half of the total area of all accreting 
cells in the disc carry material into a single hotspot, then 
half of the mass accretion rate is carried from the disc to this 
hotspot. In such a way the accretion rate into each hotspot is 
different and depends on the structure of the field connecting 
the star to the disc. Simulations suggest that both the mass 
and energy flux will be greatest at the centre of the hotspot 
( Romanova et al.|2004 l. 

In order to calculate the variation in the density struc- 
ture of the accretion columns of V2129 Oph we assume that 
the field is not significantly distorted by the flow and that 
the accreting gas flows along the path of the field (v is par- 
allel to B) so that 
pv 

— — const. (3) 
B 

where p is the mass density (see e.g. |Mestel|[T96 8). If an 
individual accreting grid cell carries material onto an area 
of the stellar surface of A* with velocity w* and density p*, 
then the mass accretion rate into A t is, 



M 



p*v t A*. (4) 

It then follows from (J3j) and Q that the density profile along 
the path of the field line may be written as, 



p(r) = 



B(r) 



M 



(5) 



B(R*) A,v(r)' 

The density profiles therefore do not depend on the absolute 
field strength, but instead on how the field strength varies 
with height above the star. The density profiles are typi- 
cally steeper than those derived for accretion flows along 
dipolar field lines since the strength of a higher order field 
drops faster with height above the star (see Fig. |4|. Gas is 
assumed to free-fall under gravity along the path of the field 
after leaving the disc with a small initial velocity of 10 kms -1 
( Muzerolle, Calvet fc Hartmann||2001 1. The range of accre- 
tion velocities with which flows arrive at the stellar surface 
is shown in Fig. [5] We assume that accretion occurs over a 
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Figure 6. The distribution of the total accretion rate into the 
northern hemisphere hotspot (solid line), the southern hemi- 
sphere hotspot (dotted line) and the various low-latitude hotspots 
(dashed line). 



range of radii, as discussed by Gregory et al. (20071. This is 



equivalent to the approach taken previously by |Hartmann, 



Hewett & Calvet (19941, Muzerolle et al. (2001 1, Syming- 



ton et al. (20051, Azevedo et al. (20061 and Kurosawa et al. 
(2006) who have demonstrated that such an assumption re- 
produces observed spectral line profiles and variability. 

In order to estimate the number density of gas within 
the accretion columns we require the temperature distri- 
bution along the path of the field lines. The temperature 
variation within the accretion columns of T Tauri stars is 
uncertain. The most comprehensive model is that of |Martin] 
( 1996 1 who considers both heating and cooling processes. 



Here we elect to take the simplest approach and assume 
an isothermal accretion flow temperature of 7000 K - a rea- 
sonable estimate for the temperature of accreting material 
(see e.g. Hartmann et al.| 1994 Martin 1996). We obtain the 
number density profiles from 



n(r) = p(r)/p,mH, 



(6) 



where tuh is the mass of a hydrogen atom and u the di- 
mensionless atomic weight. Fig.|4]show sthe variation of the 
number density along the paths of a selection of accreting 
field lines, with those obtained for dipolar field lines shown 
for comparison. 

The changing nature of the accreting field lines as the 
source surface is increased becomes very apparent when we 
consider the path taken by the accretion flows. Fig. |6| 
shows the fraction of the mass accretion rate that impacts on 
the different regions of the stellar surface. When the source 
surface is very small, the accreting flow takes place along 
the complex smaller-scale field close to the stellar surface 
and impacts at a range of locations. As the source surface 
becomes progressively larger, however, the accretion takes 
place along the simpler large-scale field until most of the 
accretion is divided equally between two accretion funnels, 
one in each hemisphere (see Fig. |7f). 
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Figure 7. Accreting field lines for V2129 Oph for a source surface at 6.8R* (left) and at 3.2R* (middle). The right-hand panel shows 
an X-ray image for a source surface at 3.2R*. In all cases the star is viewed at rotation phase 0.8 where the maximum in the accretion 
signatures is seen. 



5 THE X-RAY BRIGHT FIELD 

For the field lines that are closed and not interacting with 
the disk we can calculate the pressure structure of the 
corona, assuming it to be isothermal and in hydrostatic equi- 
librium. Hence for a stellar rotation rate ui, the pressure at 

any point is p = poe kT J 3sds where g s = (g.B)/jBj is the 
component of gravity (allowing for rotation) along the field 
and 



9(r,0) 



-GM+/r 2 + ui 2 r sin 2 6, uj 2 r sin 9 cos 



(7) 



At the loop footpoints we scale the plasma pressure po to the 
magnetic pressure such that po{8, <f>) = KBq(8, <j>) where K 
is a constant that is the same on every field line. By scaling 
K up or down we can scale the overall level of the coronal 
gas pressure and hence the density and emission measure. 
The value of K is a free parameter. Following Jardinc et al. 
( 2006 1 we choose the value of K that gives the best fit when 



used to model the coronal emission of the stars in the COUP 
dataset. The plasma pressure within any volume element 
of the corona is then set to zero if the field line through 
that volume element is open, or if at any point along the 
field line the plasma pressure is greater than the magnetic 
pressure i.e. > 1. From the pressure, we calculate the 
density assuming an ideal gas and determine the morphology 
of the optically thin X-ray emission by integrating along 
lines of sight through the corona. 

This technique has been used for some time now to ex- 
trapolate the coronal field of the Sun. Comparisons of the 
resulting morphology of the corona with LASCO CI im- 
ages show good agreement for much of the corona, except 
near the polar hole boundaries where the potential field ap- 
proximation is likely to break down at the interface between 
open and closed field regions ( [Wang et al.|1997[ ). |Wang et al 



| |1997p found the best agreement with observations was ob- 
tained with a scaling of po oc -Bq' 9 . The question of the 
optimal scaling is one that is often addressed in the context 
of the heating of the solar corona (see |Aschwanden| ( |2001[ ) 
and references therein) and indeed recent work suggests that 
scaling laws developed on the basis of detailed solar obser- 
vations may be extrapolated to more rapidly rotating stars 



on surface magnetograms obtained using Zeeman-Doppler 
imaging has been developed extensively for rapidly-rotating 
main-sequence stars (Jardine et al. 2002a b Mclvor et al. 
2006[ ). Comparison of the predicted coronal X-ray emission 
with Chandra line shifts and emission measures confirms the 
usefulness of this technique ( Hussain et al.|[2007 1 . 

One significant difference, however, between these main- 
sequence stars and T Tauri stars is the presence of active 
accretion. The disk may truncate the corona, forcing field 
lines that might otherwise have been closed (and hence X- 
ray bright) to be open or accreting, and hence X-ray dark 
(Jardine et al. 20061. In addition, the accretion funnels 



themselves may scatter some of the X-ray photons from the 
corona, causing not only a drop in the magnitude of the X- 
ray luminosity, but also an increase in its rotational modu- 
lation. Using realistically-complex magnetic field structures, 



Gregory et al. ( 2006 2007 1 have shown that for the actively- 



( Schrijver & Title 2001 1. The use of field extrapolation based 



accreting T Tauris in the COUP sample, the observed drop 
in L x (and its rotational modulation) could be explained by 
the presence of dense accretion funnels. 

In order to test this hypothesis with the derived mag- 
netic structure for V2129 Oph, we have calculated the den- 
sity profiles along the accreting field lines and determined 
the X-ray emission measure as a function of rotational phase 
for a range of models with different source surfaces. We chose 
a coronal temperature of 3xl0 7 K which is typical for stars 
in the COUP database. A lower temperature would have 
given a lower emission measure, but the behaviour as the 
source surface is changed would have been qualitatively un- 
changed. The values of the emission measure, while high for 
most solar-mass T Tauri stars, are typical for a star of this 
mass. The average value for the case of a source surface at 
6.8R* is 7.6 x 10 53 cm 
minosity of 6.7 x 10 30 ergs 
published value of 2.5 x 10 30 ergs -1 which was derived from 
ROSAT observations which were sensitive to rather cooler 
plasma ( jCasanova et al.|1995[ |. 

The effect of the accretion funnels can just be seen in the 
X-ray image in Fig. |7| which is viewed at rotation phase 0.8 
which coincides with the inferred location of the accretion 
hotspots. The effect on the X-ray emission is rather small 



which translates to an X-ray lu- 
1 . This compares well with the 
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Figure 8. Calculated emission measure for various source sur- 
faces when scattering by accretion funnels is taken into account 
(top) or neglected (bottom). The shading demonstrates the range 
of emission measures over one stellar rotation. 



because of the relatively low density of the accretion funnels, 
and the fact that the magnetic geometry of V2129 Oph is 
such that there is a ring of X-ray bright field at high lati- 
tude which is barely obscured by the channelled funnel of 
accreting gas. Fig. Q shows the variation in the density for 
some sample field lines, for a model with a source surface 
at 6.8R* and one at 2.5R*. In each case, the corresponding 
density profile for dipolar field lines is also shown. The larger 
the source surface becomes, the more dipolar do the density 
profiles become. In both cases, however, the density close 
to the stellar surface is similar and is of a magnitude such 
that scattering from the accretion funnels can reduce the 
X-ray luminosity, although the predicted rotational modu- 
lation for this moderate accretor is only of order 10% at 
most. Fig. Q shows the corresponding X-ray emission mea- 
sure for this set of models when the accretion funnels are 
present and also when they are neglected. In both cases, in- 
creasing the source surface initially increases the emission 
measure as the volume of the closed corona increases, but 
once the source surface extends beyond the co-rotation ra- 
dius, the maximum extent of the closed corona is limited to 
the co-rotation radius and so the emission measure becomes 
independent of the source surface. 




53 7 I i i i I i i i I i i i I i i i I i i i 
0.2 0.4 0.6 0.8 1 
rotation phase 



Figure 9. Rotational modulation of the emission measure. The 
dotted lines are for a source surface at 6.8-R+ and the solid lines 
are for a source surface at 2.5-R*. Black lines show the case where 
scattering by accretion funnels is neglected, while blue lines show 
the case where scattering by accretion funnels is considered. 

The degree of rotational modulation of the X-ray emis- 
sion measure also varies with the value of the source surface. 
As Fig. Q shows, however, the octupolar structure of the 
field results in a ring of emission at high latitudes. Due to 
the inclination of the star (i = 45°) this ring is almost al- 
ways in view as the star rotates, and hence, as can be seen 
from Fig. Q, the rotational modulation is small. When the 
source surface is placed at larger radii (in this case close to 
the co-rotation radius) both the magnitude and the rota- 
tional modulation of the X-ray emission are greater. In the 
case where scattering due to the accretion funnel is consid- 
ered, the location of the minimum also varies slightly, being 
closer to phase 0.8 when the source surface is larger. The 
difference is simply due to the alteration in the shape of the 
field lines as the source surface is changed. In this case, the 
field lines that are accreting connect to the star at slightly 
higher latitudes and later phases when the source surface is 
moved out. 



6 RESULTS 

From the surface magnetogram of V2129 Oph, we have ex- 
trapolated the coronal magnetic field, and from this calcu- 
lated the location, density and velocity of the accretion flows 
and also the coronal density distribution and X-ray emis- 
sion measure. In doing this we have two free parameters. 
One is the constant of proportionality K that determines 
the plasma pressure at the base of the corona, and hence 
the overall magnitude of the X-ray emission measure. In the 
absence of observations of V2129 Oph that could determine 
these parameters, we choose to select the value of K that 
gives the best fit to the range of emission measure of the 
stars in the COUP sample. The other free parameter of our 
model is the value of the source surface where the outward 
pressure of the coronal gas is assumed to have opened up all 
of the magnetic field lines. We have taken a range of values 
for the source surface, from one well inside the co-rotation 
radius (which ensures that all field lines at co-rotation are 
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open) to one well beyond the co-rotation radius (which en- 
sures that there is a mixture of open and closed field lines 
at co-rotation). 

The principal difference in the accretion process as we 
increase the value of the source surface is that accretion 
takes place along progressively simpler field structures and 
impacts on progressively fewer sites at the stellar surface. 
We find that in order to have the majority of the accretion 
impacting on a single high-latitude site in each hemisphere, 
we need to place the source surface beyond approximately 
7R«, close to the Keplerian co-rotation radius at 6.8R*. At 
this distance from the star, the magnetic field lines that 
intersect the disk have a strength of only a few Gauss and 
are certainly weak enough to be overwhelmed by the ram 
pressure of the accreting gas. 

As the source surface is increased and the extent of the 
closed-field corona increases, we find that the magnitude of 
the predicted X-ray emission measure increases, and its ro- 
tational modulation also increases. Rotational modulation 
can be caused either by self-eclipse of X-ray bright regions, 
or by scattering from the accretion funnels themselves. For 
V2129 Oph, this predicted rotational modulation is small 
(only around 10% at most) partly because the X-ray bright 
regions remain in view as the star rotates, and partly be- 
cause it is only a moderate accretor and so the density in 
the accretion funnels is relatively small. 



7 CONCLUSIONS 

V2129 Oph is a young star in which a small radiative core 
appears to have already formed. For this star, the dynamo 
produces a field that is predominantly antisymmetric, with 
a complex surface distribution of polarities in which the 
strongest mode is that of a 1.2kG octupole. A weaker dipole 
component is also present, with polar strength of 0.35kG . 
By extrapolating this surface field into the corona, we have 
been able to determine the form of the field at a range of 
distances from the surface. The observation that accretion 
takes place mainly at one range of longitudes, into a re- 
gion of positive polarity at the surface, suggest that at the 
field at the inner edge of the disk is much simpler than the 
surface field. This constrains the rate at which the field com- 
plexity decays with height and sets the limit of the closed- 
field corona at a height close to the co-rotation radius of 
6.8R* which is consistent with traditional disk-locking mod- 
els for star-disk interaction (Konigl 1991 Collier Cameron & 
Campbell|1993 \ . This is also consistent with the predictions 



of Jardine et al. ( 2006 1 that for a T Tauri star of this mass, 



the closed-field corona should extend out to approximately 
the co-rotation radius. The fact that such a relatively weak 
dipolar component is able to trigger accretion as far out 
as the co-rotation radius remains a challenge for accretion 
models which have traditionally considered purely dipolar 
fields. Progress is being made, however, in modelling more 
complex field topologies ( Long et al.|2007 1. 

Thus, although Fig. [7|shows that the X-ray emission is 
concentrated within a height of 1R* above the surface, the 
closed-field corona extends out much further. The degree 
of field complexity that might be inferred by tracing the 
X-ray bright loops is much greater than that seen in the 
accreting loops. This has parallels in the solar case, where 



the X-ray corona is compact and very highly structured, 
while the K-corona, which is also composed of gas at 10 6 K 
trapped on closed field lines, is much smoother and extends 
out to typically 2.5Rq. The solar K-corona is not detected 
in X-rays because the density of the gas at these distance 
is too low. In a similar manner, V2129 Oph also appears to 
produce a large-scale, simple corona which is not apparent in 
the X-ray images, but which is in some sense "illuminated" 
by the accretion flows along it, without which it may have 
passed undetected. 

More observations of T Tauri stars of different masses 
and evolutionary states are needed to determine if this field 
structure is particular to young stars with radiative cores, or 
if it is also seen on fully-convective stars which may produce 
magnetic fields through a different type of dynamo. Simul- 
taneous phase resolved X-ray and spectropolarimetric ob- 
servations are also needed to determine if X-ray light curve 
eclipses are indeed visible in accreting systems and if they 
occur at the phases predicted on the basis of the derived 
magnetic field structures. 



REFERENCES 

Altschuler M. D., Newkirk, Jr. C, 1969, Solar Phys., 9, 131 

Aschwanden M., 2001, ApJ, 560, 1035 

Azevedo R., Calvet N., Hartmann L., Folha D. F. M., 

Gameiro F., Muzerolle J., 2006, A&A, 456, 225 
Beristain G., Edwards S., Kwan J., 2001, ApJ, 551, 1037 
Bessolaz N., Zanni C, Ferreira J., Keppens R., Bouvier J., 

2008, A&A, 478, 155 
Bouvier J., Alencar S. H. P., Harries T. J., Johns-Krull 

C. M., Romanova M. M., 2006, ArXiv Astrophysics e- 

prints 

Casanova S., Montmerle T., Feigelson E. D., Andre P., 

1995, ApJ, 439, 752 
Collier Cameron A., Campbell C. G., 1993, A&A, 274, 309 
Daou A. G., Johns-Krull CM., Valenti J. A., 2006, A J, 

131, 520 

Donati J.-F., Jardine M. M., Gregory S. G., Petit P., Bou- 
vier J., Dougados C, Menard F., Cameron A. C, Harries 
T. J., Jeffers S. V., Paletou F., 2007, MNRAS, pp 760-+ 

Eisner J. A., Hillenbrand L. A., Carpenter J. M., Wolf S., 

2005, ApJ, 635, 396 

Favata F., Flaccomio E., Reale F., Micela G., Sciortino S., 
Shang H., Stassun K. G., Feigelson E. D., 2005, ApJs, 160, 
469 

Flaccomio E., Micela G., Sciortino S., Feigelson E. D., 
Herbst W., Favata F., Harnden F. R., Vrtilek S. D., 2005, 
ApJs, 160, 450 

Gregory S. C, Jardine M., Cameron A. C, Donati J.-F., 

2006, MNRAS, 373, 827 

Gregory S. G., Wood K., Jardine M., 2007, MNRAS, 379, 
L35 

Giidel M., Schmitt J., Benz A., Elias II N., 1995, A&A, 
301, 201 

Hartmann L., Hewett R., Calvet N., 1994, ApJ, 426, 669 
Hussain G., Brickhouse N., Dupree A., Jardine M., van 

Ballegooijen A., Collier Cameron A., Donati J.-F., Favata 

F., 2005, ApJ, 621, 999 
Hussain G. A. J., Jardine M., Donati J.-F., Brickhouse 



10 MM. Jardine et al. 



N. S., Dunstone N. J., Wood K., Dupree A. K., Collier 
Cameron A., Favata F., 2007, MNRAS, 377, 1488 

Jardine M., Cameron A. C, Donati J.-F., Gregory S. G., 
Wood K., 2006, MNRAS, 367, 917 

Jardine M., Collier Cameron A., Donati J.-F., 2002a, MN- 
RAS, 333, 339 

Jardine M., Wood K., Collier Cameron A., Donati J.-F., 

Mackay D. H., 2002b, MNRAS, 336, 1364 
Johns-Krull C. M., 2007, ApJ, 664, 975 
Johns-Krull C. M., Valenti J. A., Saar S. H., 2004, ApJ, 

617, 1204 

Getman, K. V., Flaccomio, E., Broos, P. S., et al. 2005, 

ApJS, 160, 319 
Konigl A., 1991, ApJ, 370, L39 

Kurosawa R., Harries T. J., Symington N. H., 2006, MN- 
RAS, 370, 580 

Long M., Romanova M. M., Lovelace R. V. E., 2005, ApJ, 
634, 1214 

Long M., Romanova M. M., Lovelace R. V. E., 2007, MN- 
RAS, 374, 436 

Marino A., Micela G., Peres G., Sciortino S., 2003, A&A, 

407, L63 
Martin S. C, 1996, ApJ, 470, 537 
Matt S., Pudritz R. E., 2008a, ArXiv e-prints, 801 
Matt S., Pudritz R. E., 2008b, ArXiv e-prints, 801 
Mclvor T., Jardine M., Mackay D., Holzwarth V., 2006, 

MNRAS, 367, 592 
Mestel L., 1968, MNRAS, 138, 359 

Mohanty S., Jayawardhana R., Basri G., 2005, ApJ, 626, 
498 

Muzerolle J., Calvet N., Hartmann L., 2001, ApJ, 550, 944 
Rebull L. M., Stauffer J. R., Ramirez S. V., Flaccomio E., 

Sciortino S., Micela G., Strom S. E., Wolff S. C, 2006, 

AJ, 131, 2934 

Riley P., Linker J. A., Mikic Z., Lionello R., Ledvina S. A., 
Luhmann J. G., 2006, ApJ, 653, 1510 

Romanova M. M., Long M., Kulkarni A. K., Kurosawa 
R., Ustyugova G. V., Koldoba A. K., Lovelace R. V. E., 
2007, in Bouvier J., Appcnzeller I., eds, IAU Symposium 
Vol. 243 of IAU Symposium, MHD simulations of disk- 
star interaction, pp 277-290 

Romanova M. M., Ustyugova G. V., Koldoba A. V., 
Lovelace R. V. E., 2002, ApJ, 578, 420 

Romanova M. M., Ustyugova G. V., Koldoba A. V., 
Lovelace R. V. E., 2004, ApJ, 610, 920 

Romanova M. M., Ustyugova G. V., Koldoba A. V., Wick 
J. V., Lovelace R. V. E., 2003, ApJ, 595, 1009 

Safier P. N., 1998, ApJ, 494, 336 

Schrijver C. J., Title A. M., 2001, ApJ, 551, 1099 

Shevchenko V. S., Herbst W., 1998, AJ, 116, 1419 

Shu F., Najita J., Ostriker E., Wilkin F., Ruden S., Lizano 
S., 1994, ApJ, 429, 781 

Siess L., Dufour E., Forestini M., 2000, aanda, 358, 593 

Stassun K. G., van den Berg M., Feigelson E., 2007, ApJ, 
660, 704 

Stassun K. G., van den Berg M., Feigelson E., Flaccomio 

E., 2006, ApJ, 649, 914 
Symington N. H., Harries T. J., Kurosawa R., Nay lor T., 

2005, MNRAS, 358, 977 
Valenti J. A., Johns-Krull C. M., 2004, apss, 292, 619 
van Ballegooijen A., Cartledge N., Priest E., 1998, ApJ, 

501, 866 



von Rekowski B., Brandenburg A., 2004, aap, 420, 17 
von Rekowski B., Piskunov N., 2006, Astronomische 

Nachrichten, 327, 340 
Wang Y.-M., Sheeley Jnr N., Hawley S., Kraemer S., 

Brueckner C, Howard R., Korendyke C, Michels D., 

Moulton N., Socker D., 1997, ApJ, 485, 419 
Yang H., Johns-Krull C. M., Valenti J. A., 2005, ApJ, 635, 

466 

Yang H., Johns-Krull C. M., Valenti J. A., 2007, AJ, 133, 
73 



